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Abstract. We calculate the three loop contribution to the /3-function of the gauge coupling 
constant in a general, anomaly-free, renormalisable gauge field theory involving a single gauge 
coupling using the background field method in the MS scheme. 
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The renormalisation group functions play an important role in governing the fundamental 
properties of quantum field theories. For example asymptotic freedom is a consequence of 
the sign of the one loop contribution to the gauge /3- function, f3g^ for a non-abelian gauge 
theory Q. The variations with scale of the coupling constants of a given theory are established 
from knowledge of the /?-functions, and so to relate physics at different scales in a given theory 
it is important to know the /3-functions as accurately as possible. For QCD, there has been 
impressive progress recently with the provision of the four loop term for fig in the MS scheme, [0]. 
This extended the scheme independent two loop results of Q and the three loop calculation of y , 
^. Evidently QCD does not, however, represent the most general renormalisable anomaly-free 
gauge field theory which one could consider in four dimensions; such a theory would involve both 
fermions and scalars, with Yukawa and quartic scalar interactions. Clearly this general gauge 
theory would contain the standard model (SM) as a particular case. The one loop /3-functions 
for this most general theory, were first written down in whilst the two loop corrections were 
calculated in the MS scheme in ^. 

Although the result of ^ represents substantial progress in refining the running of the strong 
coupling constant, it transpires that the renormalisation of the most general gauge theory has 
not yet been performed at three loops. In this letter, therefore, we extend the result of [Q, ^ for 
I3g to include Yukawa and quartic scalar interactions. We perform the calculation in such a way 
that the expressions for both the conventional dimensional regularisation MS (or MS) scheme 
and (in the super symmetric case) the dimensional reduction, (DRED), scheme |9[ can be easily 
extracted. In the latter case we can compare our result with that of This is especially 

worthwhile since the result of was arrived at indirectly in the non-abelian case. Therefore, 
in fact, our results will include the first explicit calculation of /Jg^'*^^^^. We defer to a later 
publication |11], the full details of the calculation where we will also provide the extension to 



the multi-gauge coupling case, the three loop MS expressions of the remaining renormalisation 
group functions of the other couplings and the results in the special case of the SM. 

Another motivation for undertaking this three loop renormalisation rests in recent exciting 
developments for constructing four dimensional conformal field theories, which are not neces- 
sarily supersymmetric, based on ideas from the Maldacena conjecture [^ ]. In [|l3|, Kachru and 
Silverstein constructed four dimensional field theories on orbifolds so that (at least in the large 
Nc limit) their /3-functions vanished, leading to (approximate) conformal symmetry. Inspired 
by this, Lawrence et al |14] conjectured that the orbifold approach could be refined so that 



the field theories were conformal for finite Nc with and without supersymmetry. Moreover, it 
was proposed in [15| that this framework could provide a non-super symmetric resolution of the 



gauge hierarchy problem. In essence, above a certain energy scale of the order of 1 — 10 TeV, 
a conformal field theory exists on an orbifold. The low energy theory (obtained by breaking 
the conformal invariance softly, by the addition of operators of dimension less than four) can 



be phenomenologically viable [^]. Subsequently ||15[-[21|, the abelian and non-abelian orbifold 
models satisfying these criteria were constructed. The method used was to impose that all 
the renormalisation group functions of the most general four dimensional gauge theory vanish, 

(2) 

thereby ensuring conformality. For instance, the results for l3g and the one loop scalar and 
Yukawa /3-functions have been analysed. This is a systematic approach since it determines pre- 
cisely the (non-supersymmetric) models which are not just conformal at large Nc but also when 
Nc is finite. A large set of possible theories emerged; however it is as yet unclear as to which 
is the strongest candidate for realistic physics beyond the SM. Moreover, it may be the case 
that the higher order corrections to the renormalisation group functions could restrict this set 
of such possible models further |18, ^]. Thus our provision of new results at three loops here is 
relevant. 
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We consider a general renormalisable quantum field theory defined by the following La- 
grangian, L: 



1 
4! 



'^abcd4>"'<P^4''^4''^ + gauge fixing + ghosts + mass terms 



(1) 



where Xabcd is totally symmetric, y"*-^ and its complex conjugate 1^" are symmetric on their ij 
indices and is the covariant derivative, 5^ — igA^. 

We extract divergences using dimensional regularisation, DREG, with modified minimal 
subtraction, MS, in d = 4 — 2e dimensions. The mass terms in the Lagrangian have been 
suppressed because they do not contribute to (3g in a mass-independent renormalisation scheme, 
such as MS. Gauge fixing and ghosts will be discussed below. 

Without loss of generality the scalar fields, are taken to be real. As a result they transform 
in an antisymmetric representation, S^, of the Lie algebra of the gauge group, G, where 



•~>ab 



ABC oC 



and f^^^ are the structure constants of G. 

We work with two-component fermions, ij^j, and their conjugates, '0-', so chiral theories can 
be considered. The fermions transform in a representation, R^, of the Lie algebra of G, where 



For cancellation of gauge anomalies we require 

Tr[i?^{i?^,fl^}] 



if 



0. 



ABCj^C 



(3) 



(4) 



In this paper, we restrict our attention to simple gauge groups and hence a single gauge 
coupling constant, g (of course the case of QED can easily be deduced as well). The extension 
to the case of several gauge couplings is non-trivial (unlike at one and two loops), and will be 



presented elsewhere |11] along with explicit results for the SM. 
The /3-function for the renormalised coupling, g, is defined as: 



Pa 



dg_ 



with go = fi^Zgg 



(5) 



where go is the bare coupling constant, and fi is the renormalisation scale. To find (3g, we 
must first calculate Zg. One might consider calculating, for example, the 3-point gauge-particle 
vertex. The graphs which contribute to this can be generated by attaching three external gauge 
legs in all possible ways to vacuum graphs; at three loops this gives rise to an enormous number 
of graphs. In addition one would of course have to calculate the gauge-particle renormalization 
factor, Za- Evaluating this to three loops involves another significant calculation, this time of 
the 2-point gauge self-energy. Moreover, difficulties arise when we consider performing the 3- 
point momentum integrals required to three loop order. We will be using the integration package 



Mincer [22| as implemented in the symbolic manipulation language Form |2S] and described 



in 1 24], which has the ability to handle only one independent external momentum and therefore 
we must set the momentum on all but two of the external legs to zero. For graphs with more 
than two external legs, this leads to the potential introduction of spurious infrared divergences. 
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which also show up as poles in e in DREG and are therefore difficult to distinguish from ultra- 
violet infinities. Whilst it is no doubt possible to perform the calculation this way, significant 
simplifications result from the use of the background field method pioneered by DeWitt in [^] 
and later extended by 't Hooft [|6| , DeWitt |^ , Boulware ^ and Abbott H , who performed a 
two loop pure Yang- Mills calculation in the Landau gauge (a = 0) , later generalised to a / by 



Capper and Maclean |3C]. More recently the method has been used to study the renormalisation 
of gauge theories with a non-semisimple colour group, [^]. 

The background field method consists of splitting the gauge field, into a quantum field, 
Q^, and a background field, B^, such that 

Af, = Q^, + B^. (6) 

Only is integrated over in the path integral and we a choose a gauge-fixing term for 
which is invariant under gauge transformations with respect to the ghost vertices may then 
be derived from this in the usual way. Corrections to the effective action are gauge invariant 
with respect to B^^ and as a result it can be shown that the following relations hold 

Za = Zq 

= 1 (7) 

Thus Zg, which we need to obtain may be obtained from the background gauge- field renor- 
malisation constant, Zb, which is easier to compute. Similarly we can deduce the gauge param- 
eter renormalisation constant, Zq,, from the quantum field renormalisation constant, Zq, and, as 
discussed below, this enables us to check our results by performing the alculation for arbitrary 
a. Consequently, both Zg and Za can be found from calculations which only involve evaluating 
2-point graphs. There are clearly significantly fewer of these graphs and they are more easily 
evaluated using Mincer as there is no need to nullify any external momenta. 

In MS, Zg and hence Zb may be written as a sum of poles in e. We define 

oo . 
n=l 

and it follows from Eqs. (||) and (^) that 

al = i^^'z-^h'. (9) 

is determined by Ai, the simple pole in Zb\ and at any given order in perturbation theory, 
the higher order poles can be calculated from the simple pole at lower orders. We have: 



and 



ftA„ + f(4 + yA + y.|, + 2A.|;)A„« 



where the reduced /3-functions, /? are independent of e and 

4 = /3g + e<7, /3y = /3y + eF, /3a = /3A + 2eA. (12) 
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Eq. (|Tl]) can be used to verify that the higher order poles generated in the final result for (3g are 
consistent with the lower order results thus providing a strong check on any calculation. 

The first step in the calculation of Zb is to generate all of the relevant graphs along with their 
appropriate combinatoric factors. We found it useful to write a package in Mathematica 
called Feynalyse, which, among other things, is capable of generating all of the required graphs 
starting from the Lagrangian, C. (As a check, we also generated the graphs using Qgraf ||3^ .) 
By producing the graphs in a MATHEMATICA environment we are then able to manipulate 
them using other routines in Feynalyse, which can draw and label each graph systematically, 
calculate combinatoric factors and signs due to fermion loops, assign Feynman rules to the 
vertices and propagators and use these to generate the necessary mathematical expressions for 
each graph. Once the complete set of graphs, signs and combinatoric factors has been found, the 
remaining calculation can be factorised into two independent calculations: the group theoretical 
factor and the momentum integration. The only vertex for which this causes a problem is the 



4-point gauge-particle interaction discussed in more detail in |11] 



The group theory is dealt with in Mathematica by using a set of specially written trans- 
formation rules which allow any group theory expression up to three loops to be reduced to a 
sum of linearly independent terms. It is important to construct such a set to ensure that all 
necessary cancellations between similar terms are carried out. We make use of the Lie algebra 
of the group generators, and R^, defined in Eqs. (Q) and (^), the Jacobi relation for the 
structure constants, f'^^'^ and the following relations 

j^Ajyaik ^ JiAiyakj ^ = 0, (13) 

which follow from gauge invariance. Most of the group theoretical factors we encounter up to 
three loops can be expressed in terms of the quadratic Casimirs of the relevant representations 
of the Lie algebra of G. These are defined as follows: 

TV(5^5^) = 5^^T{S) Stsi = C{SU (15) 
Tr{R^R^) = 5^^T{R) R^Rf^ = C{R),^ (16) 
jACD^BCD ^ s^'^'CiG) = r. (17) 

As mentioned above, we employ the anomaly cancellation condition, Eq. (Q). This justifies our 
use of the naive 75 matrix in defining two component spinors, because, at three loops, graphs 
containing two one loop fermion triangles are the only graphs which can generate non-trivial 
corrections to the gauge propagator involving the e'^'^^'^ tensor. As explained in detail in [pT]], 
Eq. (^) means that we can set e'^'^'"^ to zero, thereby avoiding inconsistencies with the 7-algebra 
in d-dimensions. 



The momentum integrations are performed using MiNCER |22, 24 1. We use it to calculate 
the divergent part of each diagram, expressed as poles in e. Although we need only calculate 
the simple pole to deduce Pg, we keep higher order poles as a check on our results. By retaining 
all the poles and the appropriate finite terms for each diagram at one, two and three loops, 
we can perform the subtraction of sub divergences in one step, at the end of the calculation, as 
explained below. First, though, we must calculate the unsubtracted divergences for all graphs 
contributing to the gauge propagator up to three loops. In order to do this using Mincer it 
is essential to label the momenta in the graph according to its 'topology', exactly as defined in 



the Mincer documentation,! 22]. Each 'topology' corresponds to a particular integration routine 
which assumes that the momenta are already correctly labelled when it is called. Due to the large 
number of graphs it is necessary to construct a program which can work out which integration 
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package needs to be called to integrate a given graph and label its momenta accordingly. The 
program we have written uses the definition of the graph in Mathematica, produced by our 
package Feynalyse, and the corresponding Feynman rules to produce a complete Form file 
containing the necessary scalar integrals along with all header files and definitions needed by 
Form and Mincer, as well as the correct call to the appropriate integration package for 
that topology. The file is completely self contained and can be read by Form which performs 
the integrations using Mincer and stores the results on disk. These are then collected and 
read into Mathematica by a further set of routines in Feynalyse, which combine the results 
with the simplified group theory factors and produce a result for the total contribution to the 
effective action from all graphs. By following this procedure we can evaluate the unsubtr acted 
contributions to the gauge propagator at one, two and three loops as a function of the bare 
couplings. However, as we are calculating a three loop result we must ensure that the Lagrangian 
has been renormalised up to two loops. This means that there are additional one and two loop 
counterterms which, when inserted in graphs with one and two loop topologies, give rise to 
contributions at the three loop level. By multiplicative renormalisability, the counterterms have 
the same form as the tree level interactions and to calculate these contributions independently 
would be extremely inefficient. Therefore, we adopt the procedure used in for the subtraction 
of sub divergences, which is summarised below. 

We begin by introducing a loop-counting parameter h, by replacing 

50 hgo, Yq hYo, Xq h? \q. (18) 

The corrections to the 2-point gauge-field interaction must be transverse, so that: 

r^. = rr(5..-^), (19) 

where q is the external momentum, and 

= ZbT^^. (20) 

Here F^^ is finite (as e — > 0) when expressed in terms of the renormalised couplings. We have 
calculated 

F^^ = l + Aoh^ + Boh^ + Coh^ + ... (21) 

where Aq^Bq,Cq - ■ ■ are functions of the bare couplings qq^ Yq, Aq, oiq. First we expand Zb in 
loops so that 

OO A 

Zb = 1 + ^ ^ = l + ah^ + hh^ + ch^ + ... (22) 

n=l ^ 

Then we re-express F^^ in terms of renormalised couplings using Eq. (P) and similarly for the 
other bare couplings. (In fact we do not need to renormalise Xoahcd as it appears for the first 
time at three loops and hence we can replace it with the renormalised value, Xabcd- We will 
however require Za at two loops and Zy, for the Yukawa couplings, at one loop.) Note that 
from Eq. (p^), a, 6, c are functions of the renormalised couplings, Y , A, a and contain only 
pole terms in e. We thus obtain 

F^^ = 1 + Ah'^ + Bh^ + Ch^ + ... (23) 



where A, B,C . . . are functions of the renormalised couplings. Now substituting Eqs. (|22|), (^ 



in Eq. (20), we simply impose that F^^ is finite as e ^ to deduce Zb. We thus have: 

a = -Pole[^] 

b = -Pole[B + aA] (24) 
c = -Pole[C + Ba + Ab], 
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where Pole [x] denotes the pole part of x with respect to e. We stress that it is important to 
keep some of the finite corrections to the lower order one and two loop graphs (ie A and B) in 
order to calculate a, b and c correctly. 

We define the loop expansion of the /3-function to be 

4 = pWh^ + pi^)h^ + p(^)h^ + ... (25) 

Similarly, Eq. ( p2|) implies that 

Ai = aih'^ + bih"^ + cih^ + . . . (26) 

where we would expect that 

ai , bi b2 ci C2 C3 

a = — , b= and c= H + (27) 

e e e"' e e'^ 

but it follows from Eq. (^) that 62 = cs = 0. Using Eq. (|10|), we find that 

/3« = -gai (28) 
= -2561 (29) 
= -3gci (30) 

Our calculations verify that 62 = C3 = and we find 

2^ 2g^T{R) g-^TjS) 
(167r )ai = (31) 

fl6 2)2. ^ 17C(G)^g^ 5C(G)g^T(fi) C(G)ff^r(g) 
^ ^ ^ 3 3 6 

g^TTjCjSf) ^ g^TT{C{R)Y-Y-) _ g^TrjCiRf) 
r It r 



aevrYci = 2857C(G)^g6 \^\^C{Gf g'^ T{R) ^ 79C{G)g^T{R) 



2 



162 162 162 

MbC{Gf g^T{S) 29C{G) g^T{R)T{S) _ C{G) g^T{Sf 

648 ^ 162 324 

1129 C(G) g^TT{C{Sf) 25 g^ TiR)Tr{G{Sf) 

216 r ^ 547 ^ 

A9g^T{S)Tr{G{Sf) 3 C(G) 5^ C(5)„b Tr(y"yf') 



216 r 4r 
7g^ C(5)2fe TV(y"y^) 2C{G) g'^ Tr(C(i?)y°y") 

6r r 
8 g4 C7(5)„fe Tv{G{R)Y''Y^) ^ 5 g^ TyjCiRfY^Y^) ^ 

+ 



3r 12 r 

g^Tr(c(i?)y"C(i?)y") _ g2 (g(g)^^ Tr(y^y"y°y^) 

12 r 6 r 

g2 C7(5)„i, Tr(y'=y'^y'=y*) g^ C7(5)^j, xr(y^y'^y'=y 



+ 



6 r 6 r 

g2 Tr(C(i?) y^y'y^y^) g^ Tr(C(i?) y«yfey<^ y^) 



24r 2r 

52 Tr(i?^y'^y''i?^y''y") 7 Tr(y'^y^)Tr(y^y"C(i?)) 

3r 12 r 



7 



12r 8r 
205 C7(G) TY{C{Rf) ^ II T{R)Tr{C{Rf) ^ 



54 r 27 r 



23 / r(5) Tr{C{Ry) ^ g'' Ti[C{RY) 29 g^ Tr(C(5)^^ 



lOSr 3r 12 r 

Agfccd ('S''^'S'^)afc {S^S^)cd g^ C{S)ah\acde\cde /„„n 

6r + 72 r ^^"^^ 



where the quartic scalar couphngs makes its first appearance in the final line |34|. We also find 



2 3 ^ mC{Gfg^ 89C{Gfg^T{R) lQC{G)g^T{Rf 
^ ^ ' ^"^ 27 27 27 

UC{Gf g^T{S) _ 7C{G) g^T{R)T{S) _ G{G) g^TjSf 

27 54 108 ^ 

IIC(G) g^TijGjSf) _ 2g^T{R)Tr{G{Sf) _ 

9r 9r 
5^ Tr{C{Rf Y"- Y"") c/^ Tv{C(R) Y" G{R) Y'') 



2r 2r 



+ 



c/2 Tr(c(i?) y^ y*) g'^ tt{g{r) y" y^ y^ y^) 

12r 3r 

c/2 Tr(c(i?) y'^ y'' y^ y") r(5) TV(C7(i?)2) 

127 187 ^ 

g2 Tr(y'^y^') TriY^Y^CiR)) g^ T{S) Tr(C(5)2^ 



6r 18 r 

\2\ o „6^^n\ rTi„/^^n\2\ 



+ 



llC(G)g6Tr(C(i?)^) _ 2g^T{R)TriG{Ry, 

9r 9r ' ^ ^ 

In any calculation it is important to check the internal consistency of the computations as well 
as ensure that they agree with any results in the literature. Given that our calculation involves 
a large degree of automatic computation it is vital that we apply all possible checks to Eqs. p^ ) 
and (|3^). First we consider internal consistency. We have verified that the gauge self-energy 
is transverse, in other words that the sum of all the relevant Feynman diagrams yields zero 
when contracted with the external propagator momentum. We have also checked that Pg is 
independent of the gauge parameter, a. Notice also that there are no C(3)-terms in Eq. (|3^); 
this non-trivial cancellation was to be expected, occurring as it does in the special case of QCD. 
Finally, we found &2 = C3 = 0, in accordance with Eq. ( |Tl| ) and, by calculating /3y to one loop 
we have explicitly verified that C2 is predicted to be exactly as given in Eq. ( p^ ) above. All this 
is evidence that we have performed the subtraction of sub divergences correctly. 

In addition to these internal consistency checks there are also several existing calculations 
of the /3-function up to and including some three loop calculations for specific cases. Many of 
these are concerned with supersymmetric theories, using DRED, and so we cannot obtain them 
easily from the DREG results given above, as this regularisation scheme violates supersymmetry. 
However, we in fact performed the computations [0], in such a manner that both DREG and 
DRED results are easily extracted. The 7-algebra was done in d = D — 2Ge dimensions and 
the momentum integrals in d = 4 — 2e dimensions. Lack of space prevents us from reproducing 
the general result here. The DREG results (which we give) are obtained by choosing D = A 
and G = 1. Choosing instead Z? = 4 and C = we can obtain the DRED result for a general 
= 1 supersymmetric theory in four dimensions. (A full justification of this is provided in ||Tl[| .) 
We need not reproduce this result here as it agrees entirely with that of p^]. As mentioned 
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earlier, in |jT^ the non-abelian result was actually deduced from the explicit calculation for the 

(3) 

abelian case. Therefore we have performed the first direct calculation of f3g for a general = 1 
supersymmetric non-abelian theory. Evidently this supersymmetric result includes the special 
cases of = 2 and = 4 theories; but as an additional check we can obtain these by setting 
C = 0, with D = 6 and D = 10 respectively, thereby reproducing the results of for these 
cases. (Of course in both cases Pg vanishes beyond one loop). Moreover, for the = 4 theory, 
but now using DREG, we obtain exactly the same non-zero result first obtained in ||3^. Apart 
from supersymmetric results, we have also verified that both the one and two loop /9-functions 
agree with those already given in the literature ||7|, ^. In addition, by setting the Yukawa and 
scalar couplings to zero, we have checked that our results are in exact agreement with the 
calculation of the QCD MS /3-function at three loops given in @, |5|, but not calculated using 
the background field method. Therefore our results pass a stringent series of tests and are in 
agreement with all existing results in the literature. 

Our main result, Eq. (|33|), involves a lot of terms compared with the one and two loop 
expressions, Eqs. (^Tj) and (p^). This is, of course, dictated by the structure of the Feynman 
diagram topologies at three loops. For instance, one can have two distinct fermion loops giving 
rise to terms with two traces over the Yukawa indices. We have expressed the results as far 
as possible in terms of the usual colour group Casimirs, C{S), C{R) and C{G); but at three 
loops two terms remain in which it is not possible to combine two generators, R^, S^, using 
their group properties, with one such term depending on the Yukawa couplings. It is interesting 
that this does not occur in the supersymmetric case. This might be a complication for the 
classification of non-supersymmetric conformal field theories jl5||-pl|. 
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